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Abstract 

We consider a simple one dimensional quantum system consisting of a heavy and a light particle 
interacting via a point interaction. The initial state is chosen to be a product state, with the 
heavy particle described by a coherent superposition of two spatially separated wave packets 
with opposite momentum and the light particle localized in the region between the two wave 
packets. 

We characterize the asymptotic dynamics of the system in the limit of small mass ratio, with 
an explicit control of the error. We derive the corresponding reduced density matrix for the 
heavy particle and explicitly compute the (partial) decoherence effect for the heavy particle 
induced by the presence of the light one for a particular set up of the parameters. 
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1. Introduction 

Decoherence has become the terminology for the irreversible suppression of interference of the 
wavefunction of a quantum system due to the interaction with an "environment" ([GJKKSZ], 
[BGJKS]). 

The usual picture of decoherence, in the simple setting of a two particle system, goes as follows. 
Suppose one has a particle M with initial wave function <p(x) = (pi(x) + (fi r (x), representing 
the superposition of two wave packets ipi, more or less supported "on the left of the origin" 
and heading to the right, and (p r supported more or less on the "right of the origin" with an 
average velocity pointing to the left. Suppose that another particle m, described initially by 
the wave packet $(y), passes by and interacts with M. Assuming a small mass ratio between 
the second and the first particle, it is conceivable that the evolution of M will not be much 
affected by the interaction, while the scattering process undergone by the particle m will depend 
strongly on the position of the heavier particle. After interaction (which is assumed to be very 
fast) one then expects that the wave function describing the state of the system is of the type 
ip(x,y) = (pi(x)$i(y) + (p r (x)$ r (y), where $/ and <3> r will have spatial supports concentrated 
in distant regions for all later times. Therefore in the configuration space of the entire system 
the entangled state will appear as the sum of two disjoint components and the possibility of 
interference of the heavy particle wave packets will be reduced. 

Notice that the reduced density matrix of the particle M has in this ideal case negligible off- 
diagonal elements. In this sense, interference has been reduced and the motion of the particle 
M has become more " classical" . That is the way decoherence plays a role in the explanation 
of the emergence of classical behaviour from quantum mechanics. 

For the relevance of the mechanism of decoherence in the classical limit of Quantum Mechanics 
in the language of Bohmian Mechanics see [ADGZ]. 

In this respect it is an interesting problem to separate "pure decoherence" from the other 
effects which an environment usually produces, which are dissipation and fluctuation. That 
is, one would like to have the motion of the system not much affected by the interaction with 
the environment, while the environment produces decoherence. It is unclear whether these 
desiderata, namely a good decoherence rate and a more or less unperturbed motion can be 
consistently fullfilled in realistic physical models ([GH]). 

Explicit models where one can rigorously establish decoherence in this sense have been worked 
out in the last years. One such model has been studied in ([DS]) where the interaction of a 
particle with the radiation field has been considered. We shall now study another interaction 
which elaborates closer the idea of scattering of light particles (the environment) off a heavy 
particle (the system) (see [JZ],[GF],[T] for similar ideas). 

We consider a very simple one dimensional model of a system (a heavy particle of mass M) plus 
environment (one light particle of mass m ) interacting via a short range force (<5— interaction). 
We consider this case as useful preparation for the treatment of a three dimensional gas of light 
particles interacting with the heavy particle, which we shall address in subsequent work. 
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We wish to stress some features that make the two-particle model with 5— interaction (which 
has the advantage of being analytically easily accessible) particularly suitable as a model for 
decoherence. 

1. there exists a simple dimensionless parameter in the problem, namely the 
fraction of the masses e = 

2. Letting e become small while keeping M fixed allows to approximate the 
solution of the Schrodinger equation of the two body problem by a scattering 
solution in which the heavy particle acts as a scattering center for the light 
one. The error is 0(e). The time scale on which this approximation holds 
is of course given by the time the light particle needs to pass the heavy 
particle. This approximation is the starting point of the analysis in [JZ]. 

3. The decoherence effect (i.e. the amount by which the off diagonal elements 
of the reduced density matrix are reduced) can be explicitely computed 
(see ( |3.2U| ) and the discussion follwing it) and is, in the relevant regime, 
of the order of aomh~~ 2 5, where 5 is the initial spread of the light particle 
and Q!o is the strength of the potential ((ctQm)~ l Ti 2 is the effective range of 
interaction) all of which can also be chosen e-dependent. 

We wish to warn the reader, that the point interaction we look at here in form of the (^-potential 
is for finite Oq not a hard core interaction. The case cto — > oo corresponds to hard core. 

The paper is organised as follows. 
In Section 2 we introduce the model and characterise the asymptotic dynamics of the two- 
particle system for small mass ratio and state the main approximation result. 
In Section 3 we show the attenuation of the off-diagonal terms in the reduced density matrix for 
the heavy particle and we compute explicitly the probability distribution for the position of the 
heavy particle, showing reduction of the interference effects with respect to the non interacting 
case. 

In Section 4 we give the proof of the main result of the paper. 

In the appendix we recall the derivation of the explicit solution of the Schroedinger equation 
of the two-body system in interaction via a delta potential in dimension one. 



2. Expression for small mass ratio 

In this section we shall study the Schroedinger equation for the two-particle system in one 
dimension described by the hamiltonian 

H = --£-A R -£-A r + a 6{r-R), a > (2.1) 

2M 2m 
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In ( |2.1|) we have denoted by R the position coordinate of the heavy particle with mass M and 
by r the position coordinate of the light particle with mass m. The interaction potential is 
chosen to be a repulsive point interaction of strenth a - 

It is well known that ( |2.1|) is a well defined positive and selfadjoint operator in L 2 (H 2 , drdR), 
which is also a solvable model ([AGH-KH]). 

In fact, for an arbitrary initial state tpQ = ipo(r,R), the solution of the Schroedinger equation 
can be explicitely written as (see [S] and the appendix) 



ij;(t, r,R) = J dr'dR'Mr', R'M (*, ~(R ~ #) + jj(r - r') 

• U$(% (r-R)- (r - R')) - ^ due-^U^t, u + \r - R\ + \r' - R'\) 

n Jo 

where we have introduced the reduced mass and the total mass of the system 

mM 



(2.2) 



v = m + M 



m + M' 

and the integral kernel of the free unitary group Uq (t) corresponding to the mass M. > 



(2.3) 



U^(t,x-x')=e-^ H o( x -x') 



' M iM, x _ xl) 2 

— r-r- e e ^ > x , x G IR 
Gtt)(U 



(2.4) 



We are interested in the case of an initial state in a product form. Then we fix two real valued 
smooth functions (for ease of formulation we assume that they are in Schwartz space S) 



(2.5) 



where || • || denotes the norm in L 2 (IR). For later use, it will be convenient to choose g even. 
Using / and g we define now the states in such a way that we can easily read of the relevant 
physical scales, i.e. we code the states by the physical parameters Rq, Pq, a, r , go, 5 as follows 



fa,Ro,P Q {R) - fa,R ,Po(fy + fa,R ,P ( R ) 



fa,Ro,P (R) 



f 



R±R \ ±i PcL 



a 



95. 



1 fr-r \ -9o r 



a, 5, R , P , qo > 0, r G IR, R > a + 5 + \r \ 



(2.6) 

(2.7) 

(2.8) 
(2.9) 
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The choice ( |2.9| ) is not essential for the most part of the paper, but it sets already a geometrical 
picture which puts the results in the right perspective (see below ( j2.10| )). Later on we shall use 
this particular choice for computing effects. Note that the spread of the wave function of M is 
not given by cr but by R . 

The initial state that we consider in the following is 



,Ro,P (R) (2-10) 

The initial state (|2.1U|) is a (pure) product state for the whole system, i.e. no correlation is 



assumed between the two particles at time zero. 

The heavy particle is assumed to be in a superposition of two spatially separated wave packets, 
one localized in R = —Ro with mean value of the momentum P Q and the other localized in 
R = Rq with mean value of the momentum — P Q . The light particle is localized around r , in 
the region between the two wave packets, with positive mean momentum go- 
To simplify the notation, in the rest of the paper we shall drop the dependence of the initial 
state on Rq, Pq, ro, go- Moreover, for the convenience of the reader, we collect here some notation 
which will be used later on 

777 f 

c = m' ^ = m M ' z/ = ( 1+e ) M ( 2 - n ) 

a = -za- (2- 12 ) 

fc = |, K=^ + k (2.13) 

T : L 2 {JR 2 ,drdR) -> L 2 (IR 2 , dxidx 2 ), 

/ M m \ 

(Th)(x 1 , x 2 ) = h(x 2 + — -xi, x 2 —J7 Xl ( 2 - 14 ) 

V m + M m+ M J 

A ± = (±Rq - a, ±R + a) (2.15) 
and finally c will denote a positive numerical constant. 

We shall now characterize the asymptotic behaviour of the wave function for small value of the 
mass ratio e for the initial state ( 2.10|) . Letting m become small, keeping M fixed, the light 



particles spreads with speeds v ~ h/ 5m and the time by which the light particle passes M is of 
the order of Rq/v thus decreases with m, so that M does not change much its position during 
the passing of m. 

The limit dynamics will hence describe a situation in which the light particle is scattered by the 
heavy one being in some fixed position, while the heavy particle moves freely. Nevertheless, we 
shall find that the free motion of the heavy particle is modified by the scattering event. In this 
heuristic argument we kept all the other physical parameters fixed except for the interaction 
strength. In fact, in order to keep the interaction effective on the light particle we need to scale 
«o in such a way that a^m ~ 0(1). There is of course no need to keep the other parameters 
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fixed, in fact one may well imagine 5 and Ro increasing with m, so that the kinetic energy of 
m stays finite and the spread of the M increases. We shall not discuss such choices here, but 
the estimates are detailed enough, so that other scalings can be easily discussed. This might 
become relevant in a model where the heavy particle is immersed in a gas of light particles. 
In order to formulate the main result of this section, we define the integral operator 



(W^' x °h)(k) = f dxh(x) (e- ikx + TZ^e-^e^-^ 1 ) , 7 > 0, x G 1R (2.16) 
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7 — i\k\ 

where the integral kernel in ( |2.16| ) is the generalized eigenfunction of the hamiltonian 

1 



(2.17) 



H. 



(2.18) 



and TZ^(k) is the corresponding reflection coefficient (see e.g. [AGH-KH]). Moreover we intro- 
duce the wave operator VL+ X ° associated to H lyX0 , explicitely given by 



(2.19) 



where h denotes the Fourier transform of h. 

With the above notation the asymptotic wave function, which will be denoted by 4> a (t), is 
explicitely characterized in the following theorem. 



Theorem 1. Given the initial state ( \2.10j ), then for any t > the following estimate holds 



\\m-nt)\\<[j + B}e 



(2.20) 



where 



r(t, r, R) = J dyf a (y)U M (t, R-y) (W^g s ) 



m 



/ dyf a (y)U™(t,R-y) (n 



9s 



(mr 



f2.21^ 



and A,B are positive, time-independent constants whose detailed dependence on the physical 
parameters characterising the interaction and the initial state will be given in section 4- 

Remark 1. Note that ip a (t, r, R) is close to what we described in the introduction. Think of f a 
as we do as consisting of two well concentrated wavepackets, then we have the light particle's 
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scattered wavefunction correlated with the two mean positions of the heavy particles, i.e. read 
in [(Q+ y )~ 1 gs] the y morally as the scattering center. 

The result of theorem 1 can be rephrased in terms of reduced density matrix for the heavy 
particle, which is defined by the integral operator p(t) in L 2 (H) given by the kernel 

p(t, R, R!) = J drip(t, r, Rj${t, r, R') (2.22) 
We also introduce the integral operator p a (t) defined by 

p a (t,R,R') = J drip a (t,r,R)^(t,r,R') 

= J dyfMKH, R ~y)j dzU(z)W{t, R' - z)i(y, A (2-23) 



where 



I(y,z) = J dk{Wl> y g 5 ){k){Wrgs){k) = ((^T^s, (n?")" 1 ^) (2-24) 



Formula (|2.24| ), obtained through heuristic considerations, has been the main ingredient in the 



description of scattering induced decoherence in [JZ]. 
Observe that, from ( ]2.23| ),( p.24|) one has 



p\t) = U^(t)p a U M (-t) (2.25) 
where Pq is defined by the integral kernel 

Po(v,z) = fMTMZM (2-26) 

It is easily seen that T(y, z) = X(z, y), \1(y, z)\ < 1 and the equality holds only if y = z. 
Then is a self-adjoint and trace-class operator, with Tr(pg) = 1; it is also positive since 



(h,p a h) = J dyh(y) J dzh(z)f a (y)f a (z) \ dk(W a «g s )(k)(W«>*g s )(k) 
dk I dyh(y)f a (y)(W a ' y g s )(k) 



(2.27) 



Moreover we have 
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Tr((p«) 2 ) = J dydz\f a (y)\ 2 \f a (z)\ 2 \l(y,z)\ 2 < 1 (2.28) 

We conclude that and its free evolution p a (t) are density matrices describing mixture states 
and by Theorem 1, for any t > 0, one has 

Tr(\p(t) ~ P a (t)\) < (j + B ) e ( 2 -29) 

This means that in our asymptotic regime the motion of the heavy particle is a free evolution. 
On the other hand the presence of the light particle has a relevant effect, since it produces a 
transition of the initial state of the heavy particle from po(y, z) = fa(y)fa(z) to p^(y, z). 
We shall see in the next section that this is the origin of the decoherence effect on the heavy 
particle. 

Finally, it is worth to mention that the dynamics of the system can be equivalently described 
by the Wigner function. From (|2.23|) , (|2.24 ) we see that the asymptotic form of the reduced 



Wigner function describing the motion of the heavy particle is the free evolution of 

W$(R, P) = ^J dxe iPx f a (R - \x)%{R + ~x)I(R -\x,R+ \x) (2.30) 



3. Size of Decoherence 



Here we discuss an application of formula (2. 23), (2. 24) to a concrete example of quantum 
evolution and we give an explicit computation of the decoherence effect. 



We shall consider the initial state (|2.10|) with the further assumptions (which are only done for 
ease of presentation) 



Co°°( 



-1,+1), 



(3.1) 



and 



a < - < Rq - |r | , 5 < Rq - |r | (3.2) 
a 

i.e. the spreading in position of the wave packets (which are superposed in a gross superposition) 
of the heavy particle is much smaller than the effective range of the interaction and this, in 
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turn, is much smaller than the separation between the two particles. Moreover the light particle 
is well separated from each wave packet of the heavy one. 
Notice that ( |3.2| ) obviously implies ^ C 1. 

Using assumptions (|3.2|) we can give an estimate of the basic object I(y, z) for y,z G and 
then we can find a more suitable expression for the reduced density matrix of the heavy particle. 
In order to formulate the result, we define the parameter 



A = J dk\g 5 (k)\ 2 T Q (k) = 1 + J dk\g s (k)\ 2 TZ a (k) 



(3.3) 



where 



ik 



1 + TZJk), 7>0 



(3.4) 



7 — ik 

is the transmission coefficient associated to a point interaction of strength 7 (see e.g. [AGH 
KH]). Then we have 



Proposition 2. Assume Then 



( 1 5 

sup \T{y,z) - 1| < c Ucr+— ri\ + ~5 Fl 

y ,z<=A± V ac(Ro -\r \) R - \r \ 



sup \2(y,z) — A\ — sup 
seA+,zeA- yeA-,zeA+ 



X{y,z)-A 



< c 



+ 



a(R - |r |) -R - ko 



(3.5) 
(3.6) 



Proof. Using the shorthand notation (3 = a5, we note that for y E A~ 



-j= {W? v g s ) (jj = e^ ko5 - k ^g(k - k S) + Up{k)e^ s+ ^^^g{\k\ + k 5) 



(3.7) 



and for y G A + 



(3.8) 



where we have used the fact that 



9s 



(Jfe) = \/~8g(k8 - ko6)e 



i(k—ko)ro 



(3.9) 
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Then for y, z G A we have 



l(y,z) =1+J dk\~g(\k\ + k 5)\ 2 \n p (k)\ 2 e-^ k ^ 
+ J dkg(k - k 5)g(\k\ + k^n^ky^^ 
+ / dkg(\k\ + k 5)g(k - k 5)n~p{k)e- l{k+ ^ )r -^ 
= 1 + r dk\g(k + k 5)\ 2 \n p (k)\ 2 (e- 2tkM ^ - l) 



+ 



o 



dkgik — k Q 5)g(k + kQ5)lZp{k)e 



+ / dkg{k + k 5)g(k - k Q 8)TZp(k)e 

j 

= 1 + di + a 2 + a 3 



-2ik- 



(3.10) 



where we have used the identity IZp + IZp + 2\Hp\ 
Using (1373) we easily estimate a\ 



and the fact that g is even. 



lad < 2 



\y- A 



dk\g(k + k Q 5)\ 2 k\TZa(k)\ 2 < 4- / dk\g(k + k 5)\ 2 k\1l3(k)\ 2 
Jo 



<2aa 

For the estimate of a 2 it is convenient to integrate by parts 



(3.11) 



a 2 



5 



< 



2|r - y\ 
_5 



00 d /— \ 
rffc— - k 5)g{k + k 8)n^k)) e 



+ / dk 

Jo 



+ dk g'(k-k 5)g(k + k 5) 
Jo 



< 



Ro - \r \ 



^-\~g(k 6)\ 2 

g{k - k 5)g(k + k S) 
g(k - k 5)g'(k + k 5) 

2 + W\\ 



Ro - \r \ Vtt 



a(R - \r \) 



(3.12) 



The term 03 is analysed exacly in the same way and then we get the estimate (|3.5| ) for y, z G A 
Since in the case y, z G A + the computation is similar we conclude that ( |3.5|) holds. 
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In order to prove ( |3.6| ) we consider the case y G A + and z G A (the case y G A and z G A~ 
can be treated exactly in the same way) and we obtain 

l(y,z) = l + J dkg(\k\ + k 5)g(\k\ - A; 5)|^(A;)| 2 e- 2i|fc|! ? +i(|fc| - fc ^ +i(|fc|+fc) 7 
+ f dkg~(k - k S)g{\k\ - k 5)Tlp{k)e- im - k)rj ^ 



+ 



+ 



J dkg(\k\ + k 5)g(k - k 8)Kp(k)e 



i{\k\+k)- 



1 + 



dk (\g(k - k Q 5)\ 2 TZ p {k) + \g(k + k Q 5)\ 2 Hp{k)) 



o 



dkg(k + k 5)g(k - k 5)\TZ l3 (k)\ e 



2-2ik- 



+ I dkg(k + k 5)g(k — k 5)TZf3(k)e 
lo 

poo 

- / dkg(k + k 5)g(k - k 5)\7lp(k)\ 2 e 



2ik- 



2-2ik- 







dkg{k + k 5)g(k — kQ5)1Zp{k)e 



-2ik- 



(3.13) 



The estimate of the last four terms of ( |3.13j ) proceeds exactly as the estimate of a 2 in ( |3.12| ). 
On the other hand 



I - j dk(\g(k- k 5)\ 2 K p (k) + \g(k + k 5)\ 2 n^(k 
dk\g(k - k 5)\ 2 * ~ %k > 



fi-ik ) 

and this concludes the proof of the proposition. □ 



(3.14) 



Proposition 2 allows us to find a further approximate form for the reduced density matrix. 



Corollary 3. Under the assumptions ( 3. and for any t > we have 



where 



Tr((p a (t)-pf(t)) 



.1/2 



< c aa + 



1 



+ 



5 



a(R - |r |) -R - \ro 



(3.15) 



pf(t) = U M (t)p f U M (-t) 



1 ^, 



A . , ™ 



(3.16) 



pl( y ,z) = -/+(*)#(*) + -fAy)fAz) + o fJ(y)fa (*) + (3.17) 
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Proof. 



Tr((p a (t)-p f (t)) 2 )=Tr((p a -pl)^ 

= \J dydz \f+(y)JJ(z)(l(y, z) - 1) + f; (y)fc(z)(l(y, z) - 1) 

+f+(y)fc(z)(I(y, z) - A) + fc(y)f*(z)(I(y, z) - A)| 2 
< sup \2(y,z) - 1| 2 + sup \l(y,z)-l\ 2 

t/,zeA+ y,zeA~ 

+ sup \I{y,z)-A\ 2 + sup \I{y,z)-A\ 2 (3.18) 
Using proposition 2 we conclude the proof. □ 

From corollary 3 and theorem 1 we conclude that the reduced density matrix for the heavy 
particle in the position representation can be approximated by the density matrix 

P f (t,R,R') = \(U A \t)f^(R)(U M (-t)JT)(R') + \{U™{t)f;){R){U™ {-tm{R') 

+^(U M (t)f^)(R)(U M (-tm(R') + ^(U M (t)f;)(R)(U M (-t)7+)(R') (3.19) 
with an explicit control of the error. 

If the interaction with the light particle is switched off, i.e. for a = 0, we have A = 1 and 
then ( |3.19|) reduces to the pure state corresponding to the coherent superposition of the free 



evolution of the two wave packets /jK 
On the other hand, if a > one easily sees that < |A| < 1 and then ( |3.19| ) is a mixed state for 
which the interference terms are reduced by the factor A and this is the typical manifestation 
of the (partial) decoherence effect induced by the light particle on the heavy one. 
The relevant parameter A (see ( |3.3|) ) is defined in terms of the probability distribution of the 
momentum of the light particle \gs{k)\ 2 and of the transmission coefficient T a (k). 
Then the decoherence effect is emphasized if the fraction of transmitted wave for the light 
particle is small. 

In particular, rescaling the integration variable in ( |3.3|) , one can also write 

A = 1 - / dz\g(z)\ 2 ,^ (3.20) 
J ad — i[z + k d) 

Then one easily sees that for ko ^> a the light particle is completely transmitted and A ~ 1, 
i.e. the decoherence effect is negligible. 
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On the other hand, for fco -C a the decoherence effect is nonzero and given by 

A ~ J dz\~g(z)\ 2 \T a5 (z)\ 2 = 1 - J dz\~g(z)\ 2 - (3.21) 

where we have used the fact that g is even. 

A further interesting question is the analysis of pf (t) in the momentum representation. 
Since momentum is a constant of the motion, the density matrix is simply given by 



h pf ° [h> n ) ~ 2h^ (J ^{n) + (lH ( n ) 

It is then clear that the decoherence effect is present also in the momentum representation and 
it is measured by the same parameter A. 

Moreover, if /+ and f~ are well separated, one easily realizes that the probability distribution 
of the momentum remains essentially unchanged with respect to the unperturbed case A = 1, 
the error being of order e. 

We analyse now the evolution in the position representation of the heavy particle exploiting 
the approximate reduced density matrix fr{t). 

We shall explicitely show that the typical interference fringes produced by the superposition 
state when the interaction with the light particle is absent, i.e. for A = 1, are in fact reduced 
when the light particle is present, i.e. for |A| < 1. 
In order to see the effect more clearly we assume 



i«4 (323 » 

The effect of the interfernce terms becomes more relevant when the supports of the two wave 
pakets Uq (t)f^ have the maximal overlapping and this approximately happens at the time 
t = r = ^pj^-- Then, from fl3.19| ), we consider 



n(r, R) = p f (r, R, R) 
1 



l(^o A V)/+)(^| 2 + l(^o^ 
Using (|3.23| ) and a standard scattering estimate (see e.g. [RS]) we obtain 



(3.24) 
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(U M (r)f+)(R) 



P (J f dxf ( x + Rp ^ e ^ x+l ^ {R ^ xf 



iuhR 



a 



P a i 
-e 



ifiRn 



\£n\\ < 



IIA/1 



2%Rq 

Proceeding analogously for (U^ (r)f~)(R) we find 



where 



n(r, R) 

ll^ilk 1 < c 



p a 



hR n 



f 



hR 



R 



/2P 

l + \A\cos(-^R + (p 



+ £i(R) 



(3.25) 
(3.26) 



(3.27) 
(3.28) 



A= |A|e 



up 



(3.29) 



For |A| < 1, formula (|3.27 ) shows that the presence of the light particle determines a reduction 
of the amplitude of the oscillations and a shift of the corresponding phases. 
Notice that the shift is negligible if ko <C a. 



4. Proof of theorem 1 

The proof of theorem 1 will be obtained through the proof of three lemmas. 
Lemma 4. Given the initial state ( 2.1C\ ), for any t > one has 

\\4>(t)-Mt)\\<Cie 

where 



(4.1) 



Mt, r,R) = J dyf a (y)UZ (t, ^R+^r-y^jJ dr'g s (r' + y)U^ (t, r - R, r') 



(4.2) 
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and 



dxx 2 / dy 



d_ 

dy 



(fa(y)gs(x + y)) 



1/2 



Proof. Using the relative and the center of mass coordinates (see ( p.!4j )), from ( |5.4j) 

(Ti/>{t))(x u x 2 ) = (UZW&WTfa) (x 1 ,x 2 ) 
where U£ Q (t) is defined in (|5.10| ) of the appendix and 



(T«( Il ,x 2 ) = /„(x 2 -^ 1 ) 9i (x 2 + ^ 1 



Moreover 



(Tifj l (t))(x 1 ,x 2 ) = J dx' 2 dx' l f a (x' 2 )g 5 (x' 1 + x' 2 )U^(t,x 2 - x' 2 )U£ o (t, x u x[) 



where 



Then we have with Ji = -f- M _ i = 

M 1+e' f 1+6 



^01 1 



||^(t)-^(t)|| 2 = ||T^(t)-T^(t)|| 2 = HTVo-TV, 112 
dx ld x 2 " " W M 

dx\dx 2 

dx\dx 2 



fa[x 2 - j^ x ij 9s [x 2 + — x iJ - U{xi)9s{x\ + x 2 ) 
In ( £ 2 - j~—; x ij #<5 ( x 2 + £i - ^T^i) ~ fa{x 2 )g 5 (x l + x 2 ) 



F ^xi, x 2 - - ^(^1^2) 

where F(xi,x 2 ) = fcr{x 2 )gs{x\ + x 2 ). By a simple Plancherel argument, we have that 



dx\dx 2 F yxi, x 2 - - — - x 2 ) 

y dari y ^ fe) ( e -^ xifc - 1 
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< dxx dk\F(x!,k)\ : 



1 + e 

from which the lemma follows. □ 



dx\%\ J dx2 



1 + e 
d 



dx' 



X\k 



F(xi,x 2 ) 



(4.9) 



Since a small value of m in the interacting unitary group U£ (i) is equivalent to a large value 
of t, in the next lemma we use a typical scattering estimate to approximate U£ (t) in (4.2). 



Lemma 5. Given the initial state ( \2.1ty ), for any t > one has 



>i(t)-Mt)\\<Y e 



(4.10) 



where 



M*,r, R) 



m 



M 



lixiht V 2niht 



g l 2ht' ^ l 2ht n 



• dr'g s (r' + e^"^' 



e ife|r-fl||r'| 

1 -i^rlr 

ao* 1 1 



(4.11) 



and 



C 2 = Cy{ /dxx 4 |/ CT (x)|* + / dx\U(x)\' 



J dyy 4 \g&(y + x )| 2 + ~^ (/ d v\9s{y 



+ x) 



1/2 



(4.12) 



Proof. We shall first estimate the difference ipi{t) — ^{t), where ^{t) is explicitely given by 



From (14.21) we have 



2iriht 

e i£ t \r-R\\r<\ 

' l-i-^:\r-R\ 



dyfMU u [t^R+^r-y 



dr'gs (r' + y) 



(4.13) 
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(Tip 1 )(t,x u x 2 ) = J dyUo(t } x 2 - y)<pi(t,xx,y) 
(pi(t,X!,x 2 ) = fa{x 2 ) dr'g s (r' + x 2 )U£(t,x 1 ,r') 



(4.14) 
(4.15) 



and analogously for ^ 2 (t) we write 



(Tip 2 (t))(x 1 ,x 2 ) = I dyUo(t,X2-y)(p 2 (t,xi,y) 
ip 2 (t,x 1 ,x 2 ) = 



(4.16) 



e^f a (x 2 ) f dr'g s (r' + x 2 ) ( e'^' - f"'^'. I ( 1.17) 



1 _ ,'M£li 

1 1 a t 



Using the isometric character of the operators T and U^t) and the explicit expression of U£ o (t) 
(see ( |5.10|) ), we have 



L (t)-i; 2 (t)\\ 2 = \\TMt)-m 



\\<fi{t)-<f 2 (t)\f 



< 2 / dx\dx 2 



f a (x 2 ) / dr'g 5 (r' + x 2 ) \ U(f(t,Xi -r') 



— Q i 2ht X l i Hi XlT ' 



2iriht 



+2 j ' dx\dx 2 f<r(x 2 ) j ' dr' 'g 5 (r' + x 2 ) f^zjr J due ^ w Uo(t, u + \xi\ + \r'\) 



2mht l - z^M 

a t 



A standard estimate for the free unitary group (see e.g. [RS]) gives 
(-0 = -T7 / dx 2 \f a (x 2 )\ 2 J dxi J dr'g & {r' + x 2 ) (e? 



= - J dx 2 \f a (x 2 )\ 2 J dxi J dr'g s (r' + x 2 ) (e 
= 2 J dx 2 \f a (x 2 )\ 2 j dr'\g s (r' + x 2 )\ 2 [ 
< \ {£) j dx 2 \Ux 2 )\ 2 J dr'r"\g s (r' + x 2 )\ 



1) e~ l -t Q 

2 



e l 2ht r — 1 



e 2 M 2 



2(l + e) 2 W 



dx 2 \f a (x 2 )\ dr'r' \g s {r' + x 2 )\ 



(4.18) 



(4.19) 
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where in the second line of ( |4.19| ) we used Plancherel theorem. 
Concerning (II), we introduce the change of variables 



jtU, Vl = ^-x x (4.20) 



and use the identity 

/ Ar>\o ant — 

• Rlasil 



dve-^—* = -4^ (4.21) 



aot 



Then 



(//) = - I dx 2 \f a (x 2 )\ 2 I dyi 



dr'gs (r' + x 2 ) e 



i\yi\\r'\ 



i-U-r' 2 



dve 



dve~ v+il ^ Mv 



< - I dx 2 \f a (x 2 )\ 2 I dyi 



e 2hta 2 Ufa ' ' — | 



+ - / dx 2 \f a {x 2 )\ 2 J dy 1 



dr'g s W' + x 2 ) e^^V^' 2 f°° dve~ v+il ^ Mv 

Jo 



dr'g s (r' + x 2 ) e i|s/l||r ' 1 (e^ 



12 1 

2ht' — I 



(///) + (IV) 



(4.22) 



The estimate of {IV) is trivial 



< /v '> £ S (&)7^'<*>''/< % + ( iL)'rf (/*' v2 l*"'' + ^)l) 2 



< 2(l + e)3 ^ t 2 / ^|/ CT (x 2 )| 2 {J dr'r' 2 \g 5 (r' + x 2 )\j (4.23) 
For the estimate of (III) it is convenient to integrate by parts the integral in the variable v 



7T 



(III) = - dx 2 \f a (x 2 )\ 2 dy 



l + 



dr'gs(r' + x 2 )e 



dve 



-u + r 
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< 



< 



2(l + e)/x 2 
nh 2 t 2 a 
2e 2 M 2 



dx 2 \f a {x 2 )\ 2 / dy 



dr'\g$(r' + x 2 )\ / dve v 
Jo 



1 + e 



-v + r 



a: 



[1 + e)Wa 



dx 2 \f a {x 2 )\' 



dr'\gs(r' + x 2 )\ 



1 + e 



a 



(4.24) 



Finally, it remains to analise the difference ip 2 (t) —ip 2 (t). Using the explicit expression of Uq (t) 
we have 



V> 2 (*,r, R)-Mt,r, R) 



m 



M 



2iriht V 2iriht 



J dr'g s (r' + le^ 



d^(0 (f^ - e^ 2 ) e -*(^fiO€ 

e ^|r-B||r'| 



a < 1 1 



(4.25) 



Exploiting Plancherel theorem and the fact that the operator ( |2.16|) is unitary one easily sees 
that 



««)-««)li 2 < 



\2HJ 



dti 4 \U0f 



(4.26) 



From (|l D ,(gl D ,(gj|) ) (gjD, (|jD , ( gjg) we conclude the proof of the lemma. □ 



In the last step we approximate ( [4.13|) using the fact that the coordinates of the heavy particle 
are slowly varying with respect to the coordinates of the light one. 



Lemma 6. Given the initial state ( \2.1(\ ), for any t > we have 

\\Mt)-r(t)\\<C 3 (t)e 

where 



(4.27) 



with 



Cs(t) 
1 



+ 



V2tt 



dzz 2 I dx 



dxdz 



iL ^ frri: ~ ' r) ' j{ ' r ' 



1/2 



a 2 + z 2 
a 2 + x 2 



dr'\az,r')\) + 



a 2 z 2 



a 2 + x 2 



C 3 (t)<^ + C 5 



dr'\r f \\((z,r')\ 



1/2 



(4.28) 



(4.29) 
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and 



i M_ f 2 



((z,r')= dtUZ)g s (r' + Z)e 



(4.30) 



Proof. From fl2.21|) and (|4.11| ) we have 



2niht V 2niht 
dr'g 5 (r' + ( e -*6<'-«V 
(V>2 - ^ a ) /r (*, r, R) + - (t, r, i?) 



c + I dr'g s {r' + i) 



e itk\r-R\\r'\ 



e "st i 



We will estimate separatly the two terms — 4> a )f r (*) an d (^2 ~ i ,a )i n W- 
Introducing the new integration variables 



m Mi? + mr 
x = — r, z = 

ht ht 



and the function f a (0 defined in (|4.30|) , we have 



a t I 



(4.31) 



(4.32) 



11 (4>2-r) fr (t)\\ 

^ dxdz 
dxdz 



(2*)* 
1 

2^ 



dU'AOe- 1 ^ dr'gs (r' + £) e 



e 



dxdz 



1 + e 



f a [z-x+ ——z )g s [x- ——z - / CT (z - x)fif(ar) 



1 + e 



I <ix<iz G ^a; — - — z, z^j — G(x, z) 



where we introduced the function G(x, z) = f a {z — x)g(x) 
Proceeding as in (|4.9|) we find 



11 - r) fr 



< 



dzz 2 / dx 



(4.33) 



(4.34) 
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For the estimate of (ip2 — i> a )in(t) we use again the change of variables (|4.32|) and we introduce 
the function ((z,r') defined in (|4.30|) 
Then we have 



< 



2-r 
i 



< 



(27T) 



dxdz 



dr'((z, t') 



e l \ x l+e 



z\\r'\ 



i\x\\r'\ 



2tt 2 
1 



2tt 2 

J2 



dxdz 



dxdz 



l-i|(l + e)z-ez 
ft 2 



1 + e)x — ez\ 1 

2 



X 



a 2 + x 2 



dr'({z, f') 



1 — - x 

a 1 

^h-Trr'Ik'l _ p i\x\\r'\ 



dr'\((z,r')\ 

2 



< 



2tt 2 



dxdz 



a 2 + z 2 
a 2 + x 2 



dr'\((z,r')\) + 



a 2 z 2 



a 2 + x 2 



dr'\r'\ \((z, r')\ 



where we have used the estimates 



1 



1 + e)x 

2 



ez\ 



\x\ 



ft 2 (\x\ - \(l + e)x 



ez 



(x — zf 



< 



ft 2 + x 2 ft 2 + ((1 + e)x — ez) 2 
2 a 2 + z 2 



and 



ft 2 + x 2 ft 2 + ((1 + e)x — ez) 2 



ft 2 + x 2 



< 



1 + e 1 



\z\ \r 



(4.35) 



(4.36) 



(4.37) 



The function ((z,r') is smooth and, using repeated integration by parts, one easily sees that 
it is rapidly decreasing when its first argument goes to infinity. The computation is long but 
straightforward and we omit the details. The conclusion is that the integral in the last line of 
( 4.35| ) is finite. Along the same line one can verify that ( 4.29 ) holds, where the constants C4, 
C5 are independent of time and then the proof of the lemma follows. □ 



Proof of theorem 1. It is an obvious consequence of lemmas 4, 5, 6. □ 



5. Appendix: explicit solution of the two-body problem 

We recall here the solution of the Schrodinger equation 

dip(t) 



ih- 



Ot 



Hi>{t), ^(0)=^o 



(5.1) 
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where H is the self-adjoint hamiltonian in L 2 (H ,drdR) given by ( j2.ll ). Using the unitary 
operator T (see (|2.14j) one obviously has 



where (rri.a^) are the relative and the center of mass coordinates 



THT- 1 = H% + H^, Hq — -—A X2 , H^ = -—A n +a 5(x 1 ) (5.2) 



Xl=T -R, * = ?^ (5.3) 

m + M 



Then the solution of (B.lj) can be written as 



= J dr'dR'Mr', 10)1% (t, ^(R-R') + ^( r - r')) U£ o (t, r - R,r' - R!) 

(5.4) 

where the interacting unitary group U£ (t) is given by 

U£ (t,x,x') = e- i: z Htl o(x,x'), x,x'eJR (5.5) 

We remark that the evolution ( |5.4j ) factorizes into a product of a free evolution in the center 
of mass coordinate and a one-body interacting evolution in the relative coordinate only if the 
initial state is of the form if)(r, R) = ^{—R + -j^r)ip2( r — R)- 

In order to compute (U^ o (t)ipo)(x) = ip(t, x) one has to solve the one-body Schrodinger equation 

ih ^l = _^_ A ( t ) + a S(x)<f(t), <p{0) = y? (5.6) 
at 2/j, 

Defining the rescaled wave function 

9(8,z) = <p(?is,^=z y j (5.7) 
one finds that 9(s) satisfies the corresponding equation with \i — h — 1 

t^ = -^A z 9(s)+a ^6(z)9(s), 9(0) = 9 , 9 (z) = ?o ( -^z) (5.8) 



The solution of ( |5.8| ) can be found in [S] 
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9(s, z) = (U£(s)6 ) (z) - a ^ J q °° dve~ ao ^ v j dz'U£(s,v + \z\ + \z'\)9 Q {z') (5.9) 

where Uq(s) is the free propagator with h = 1. Noticing that (U£ (t)ip )(x) = ip(t,x) = 
9(jr, ^x) one has 

(UZ o (t)w)(x) = (US(t)<po)(x) - ^ f due-^ u J dx'UZ(t,u+\x\ + \x'\)^ {x') 

(5.10) 



Using (|5.4j) and ( |5.10| ), we finally obtain the complete solution ( |2.2| ) of the Schroedinger equation 
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